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Abstract 

An estimate for the classical action I c iassicai f° r a Bianchi Vlh 
homogeneous spatially closed model with h = —1/9 is given by 
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where b and ao are parameters of the model, b is zero if and only if 
the relative rotation of inertial frames and matter is zero, cto is an 
initial value of an anisotropy of the expansion rate at r = tq, L* is 
the Planck length, and r m <i is a constant of integration that gives the 
maximum size of the universe for the isotropic (b = 0, ao = 0) case. 

It is assumed that the equation of state is p = (7 — l)p, where 
p is pressure and p is density. It is assumed that 7 has a constant 
value of 7 = 4/3 (to represent a relativistic early universe) for r < 77 
and a constant value of 7 = 1 (to represent a matter-dominated late 
universe) for r > 77. The approximation is valid for b small enough 

that \I c lassical[b] ~ I c lassical[b = 0] | < k . 
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An explanation for why our inertial frame seems not to rotate 
relative to the stars is found in a straightforward application of semi- 
classical approximations to quantum cosmology. Application of a sad- 
dlepoint approximation leads to the result that only those classical 
geometries whose action I c i assica i satisfies \I classical ~ haddiepoinA < ft 
contribute significantly to the integration to give the present value of 
the wave function. Using estimates for our universe implies that only 
those classical geometries for which the present relative rotation rate 
of inertial frames and matter are less than about 10 -130 radians per 
year contribute significantly to the integration. This is well below the 
limit set by experiment. The result depends on the Hubble distance 
being much larger than the Planck length, but does not depend on 
the details of the theory of quantum gravity. 



1 Introduction 



Although Newton recognized that inertial frames seem not to rotate relative 
to the stars, he seems to have taken that as evidence for the existence of 
absolute space. Ernst Mach [[Mach 1872| , [Mach 1933H was probably the first 
physicist to recognize that such an apparent coincidence requires an explana- 
tion. He based his arguments on the observation that only relative positions 
of bodies are observable. Mach further suggested that matter might cause 
inertia. 

Einstein ||Einstein 1918|| tried to include what he called "Mach's Princi- 
ple" in General Relativity, and although it is generally agreed that matter is 
a source of inertia in his theory, it is not the only source, because there are 
solutions (such as empty Minkowski space) that have inertia without matter. 
Further, solutions of Einstein's field equations include cosmologies where in- 
ertial frames rotate relative to the bulk of matter in the universe, so that an 
explanation for why our inertial frame does not rotate relative to the stars 
is still needed. 

It has been suggested (e.g. [Wheeler 1964|| ) that Mach's principle be used 
as a boundary condition to eliminate those solutions of the field equations 
that have inertia from sources other than from matter. Similar suggestions 
[ Al'tshuler 19671 , [Lynden-Bell 13571 , |Sciama et al. ~TO| , Oilman 1970 



e.g. 



[Rainc 1975|| ) that such selection might take place automatically if a sat- 
isfactory integral formulation of Einstein's field equations might be found 
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would eliminate the need for explicit boundary conditions. Other sugges- 
tions in which gravitation (including inertia) is represented by a theory 
analogous to Newtonian gravitation or Maxwell theory (e.g. [ fSciama 1953| , 



Lynden-Bell 19921 ) are also intriguing, and have the similar property that the 



gravitational field would be determined solely by the matter distribution. 

Even if we could successfully find a theory along the lines of those men- 
tioned above, in which matter somehow determined the gravitational field, 
it might not be satisfactory. Rather than explain why our inertial frame 
seems not to rotate relative to the stars, such a theory would simply impose 
that condition, and would take away all degrees of freedom from the grav- 
itational field. We are not so restrictive with the electromagnetic field, for 
example. We allow arbitrary initial values on the electromagnetic field that 
are consistent with Maxwell's equations. We don't require that the electro- 
magnetic field be completely determined by charges and currents. Just as 
the electromagnetic field should be as fundamental as charges and currents, 



the gravitational field should be as fundamental as matter [[Earman 1993 
Kuchaf 1993] , [Raine 1993 



Allowing arbitrary initial conditions for the gravitational field (consistent 
with the field equations) is inconsistent with trying to explain why our inertial 
frame seems not to rotate relative to the stars, at least on the classical level. 

On the quantum level, however, we might imagine that somehow the 
selection takes place automatically through wave interference, and this turns 
out to be the case. To show that requires calculating the action for a classical 
cosmology as a function of the appropriate parameters of the model, which 
is the goal here. 

2 An example from ordinary wave mechanics 

To help explain the ideas that follow, we first consider elementary wave 
mechanics as an example. If we have an initial single-particle state specified 
by an initial wave function < Xi,ti\tp > at time t\ then the wave function 
< ^2, £2 IV* > at time t 2 is [[Feynman and Hibbs 19651 , P- 57] 



/oo 
< x 2 , t 2 \xi,t 1 >< xi,ti\tp > dxi, (1) 
-oo 

where < a^^ki^i > is the propagator for the particle to go from (xi,ti) 
to (^2,^2)- We consider the case where the semiclassical approximation for 
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the propagator is valid. That is, fFeynman and Hibbs 1965| , p. 60] 

< x 2 ,t 2 \ Xl ,ti >^ f(t u t 2 )e^ x ^> x ^\ (2) 

where Id[x 2 , ^2! ti] is the action calculated along the classical path from 
(xi,ti) to (x 2 ,t 2 ). Thus, (|ip becomes 



< x 2 ,t 2 \i> >- f(hM / e ^' [x2 ^ Xl ' tl] < x u t^ > dx i- 



(3) 



Notice that because of the initial wave function we have an infinite number 
of classical paths contributing to each value of the final wave function. 

There are two cases to consider. In the first, I c i is not a sharply peaked 
function of X\. In that case, there will be contributions to the wave function 
at t 2 from classical paths that differ greatly from each other. 

In the second case, which we now consider, I c i is sharply peaked about 
some value of x\, say x sp . That is, we have 



_d_ 

dx\ 



Icl[X 2 , t 2 \ Xi, ti]\ Xl=x — 0. 



(4) 



Thus, X g<p IS 3j saddlepoint of the integral @, and significant contributions 
to the integral are limited to values of x\ such that 



\X\ — x 



sp I 



< 



2h 



■§jzl c i[ x 2i h] Xi, ti]\ xl=Xsp 



(5) 



If < X\, ti\ip > is nearly constant over that range, then we can take it outside 
of the integral. A saddlepoint evaluation of the integral then gives 



< x 2 ,t 2 \ip > w f(t h t 2 ) < x sp ,tt\ip > 

2-Kih 



" Id [x 2 , t 2 ] Xi, ti]\ 



1/2 



,J[Icl[x2,t2;Xsp,tl] 



(6) 



We notice from (|4]) that the momentum at ti at the saddlepoint is zero. 
That is, 

Pi\x 1= x sp = 0. (7) 
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However, for the paths that contribute significantly to the integral in (|3|), 
there is a range of momenta, namely 



\pj\ < 2h 



(J 2 



dx{ 



Id[ x 2i h'i ti] 



(8) 



consistent with (|3p and the uncertainty relation. 

As a check, using a special case, we consider the free-particle propagator 
b'cynman and Hibbs 1965, p. 42] 



< x 2 ,t 2 \x 1 ,t 1 >- 
and we choose 



m 



nl/2 



2mh{t 2 -t x ) 



exp 



im[X2 — x\ t 

2h(t 2 - h) 



< x 1} ti\ip >=< A,tx\ip > exp[—B(x 1 - A)' 



(9) 



(10) 



to represent a broad initial wave function. For this case, the integral in (j^) 
or (H) can be evaluated exactly to give 



< x 2 ,t 2 \ip >=< A,ti\i/j > 



2Bh(t 2 - h) 



im 



-1/2 



exp 



-B(x 2 - A) 2 

2Bh(t 2 -t 1 ) 



The condition that the initial wave function is slowly varying is now 



\B\ < 



in 



2h(t 2 - ti) 



so that ( PH ) is approximately 

< x 2 ,t 2 \ip >=< x 2 ,ti\ip > 
in agreement with (H), since 

m (x 2 — xi) 2 



I[x 2 , t 2 ; Xi,ti] — I c i[x 2 , t 2 ; xi, ti 



2 t 2 -h 



and 



x S p x 2 . 



(11) 
(12) 

(13) 

(14) 
(15) 
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Notice that it is the sharply peaked action that determines which classical 
paths in (H) dominate the integral in this case, not the maximum of the initial 
wave function. 

The calculations can clearly be generalized to two or three dimensions. 
The main point is that whenever the initial wave function is broad and the 
action of the classical propagator is not a sharply peaked function of x\, some 
very different classical paths may contribute to the wave function in the final 
state, even when a semiclassical approximation is valid for the propagator. 

When the classical action is sharply peaked as a function of the coor- 
dinates of the initial state, however, only a narrow range of classical paths 
contribute significantly to the wave function in the final state. This is thus a 
mechanism for selecting classical paths in wave mechanics. As we shall argue 
in the next sections, this principle has broader application. 



3 Quantum cosmology 

In the case of quantum cosmology, we have a formula analogous to (]T|) to 
give the wave function over 3-geometries g 2 and matter fields 02 on a 3- 
dimensional hyper surf ace S 2 . 

< 92, 02, S 2 \ip >= J < g 2 , 02, S 2 \gi, 0i, Si >< gi, 0i, Si|-0 > D(g 1 )D(cj) 1 ), 

(16) 

where < gi, 0i, Si\if) > is the wave function over 3-geometries g± and matter 
fields 0i on a 3-dimensional hypersurface Si, and < g 2 , 02, S 2 |gi, 0i, Si > is 
the amplitude to go from a state with 3-geometry gi and matter fields 0i on a 
surface Si to a state with 3-geometry g 2 and matter fields 02 on a surface S 2 
[ |Hawking 1979|| . D(g\) and -D(0i) are the measures on the 3-geometry and 



matter fields. The integration is over all initial 3-geometries g\ and matter 
fields 0i for which the integral is defined. 



4 Semiclassical approximation 



As in section 2, we want to consider the case where the semiclassical approx- 
imation for the propagator is valid. That is, | Gerlach 1969| 



< 22,0 2 ,S 2 |2i,0i,Si >^ /(^,02,S 2 ; (?i,0i,Si)e^^^^ lA ^, (17) 
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where the function outside of the exponential is a slowly varying function 
and I c i is the action for a classical 4-geometry. Substituting ( |P7|) into (|16|) 
gives 

< 92, 02, S 2 \^> = J f(g 2 , 2 , S 2 ; g u 1; S^iW.fe^WxA] 

< J1 >i ) S 1 |^>%)%). (18) 

Each value of the integrand in ( fL8f) corresponds to one classical 4-geometry. 
As in (|3|), there will be an infinite number of classical 4-geometries that con- 
tribute to each value of the final wave function. Here, however, we do not 
have only one single integration, but an infinite number of integrations, be- 
cause the integration is carried out over all possible 3-geometries and all 
matter fields on the initial surface. 

In the simple example in Section 2, there were two cases to consider for 
the single integration being carried out. In the first case, the classical action 
was not a sharply peaked function. In the second case, the classical action 
was a sharply peaked function so that a saddlepoint approximation could 
be applied to the integration. Following that strategy, we would need to 
consider those two cases for each of the infinite number of integrations in 



Here, however, we consider only the case where J c / is a sharply peaked 
function of g\ and matter fields 0i for each of the infinite number of integra- 



tions in (18|). We consider this case in the following section. 



5 Saddlepoint approximation for the integral 
over initial states 

We consider the case here where I c i is a sharply peaked function of g\ and 
matter fields 0i for each of the infinite number of integrations in (|TE|). In that 
case, we can formally make the saddlepoint approximation for each of the 
integrations in (113). In analogy with (11), we have the saddlepoint condition 



d 

tt- Id[92,(p2, S 2 ; #1,01, Si] 

dgi 



= (19) 

gi=gs P 
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and 

d 

nT-Ici\92, 02, S 2 ; 9i, 0i, Si] 



0, (20) 



where the derivatives in (19) and (]20f) are with respect to each parame- 
ter that defines the 3-geometry g\ and matter fields 0i. We consider the 
case where there is only one solution to the saddlepoint conditions (|19|) 
and ( ppp . In that case, (|I9D selects a single classical 4- geometry. How- 
ever, there will be a range of classical 4-geometries in the neighborhood that 
contribute significantly to the integral in (|I£|). These are determined by (e.g. 
Uourant and Hilbert 1953fl ) 



\hi[92,(p2,S2]gi,(t>i,S 1 ] - I d [g 2 ,(p2,S2]gsp,4 > sp,S 1 \\ <K. (21) 

We can formally write the saddlepoint approximation to the integration in 
©as 

< g2, 4>2, S 2 \^P > = f(g2,4>2,S 2 ;g sp ,(psp,S 1 )<g sp ,4> sp ,S 1 \ij> 

fi(92, 02, S 2] g sp , sp , SJefa^Sw^M (22 ) 



where classical 4-geometries that contribute significantly to fl2~2"|) (through 
the function f\) lie within a narrow range specified by (|2l|) . 

Equation ([H]) requires that the momentum canonical to the initial 3- 
geometry for the classical 4-geometry at the saddlepoint be zero. That is 



7T^ 



= 0. (23) 

gi=9s P 



(The extrinsic curvature on Si will therefore also be zero at the saddlepoint.) 
However, there will be a range of initial canonical momenta and a range of 
initial 3-geometries corresponding to the range of classical 4-geometries that 
satisfy (0), so that the uncertainty relations between initial 3-geometries 
and their canonical momenta are satisfied. 

Whether there is a narrow or broad range of classical 4-geometries that 
satisfy ( PT| ) depends on the second derivative of the action with respect to 
the initial 3-geometry. 

We can take the action to be 

I = J (-9 {4) ) 1/2 (L geom + L matter )d'x + i- J (gWy/ 2 Kd 3 x, (24) 
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where [ |York 1972] , |Hawking 1979|| show the importance of the surface term. 



Hawking 1979 | also points out a potential problem in that the action can be 



changed by conformal transformations, but suggests a solution. 

K = gVMKn (25) 

is the trace of the extrinsic curvature. Although the extrinsic curvature is 
zero on S\ at the saddlepoint, it will be nonzero in a region around the 
saddlepoint. The extrinsic curvature is given by 

K« - - 5 -gL, (26) 

(3) 

where g\j is the 3-metric. In this example, we take the Lagrangian for the 
geometry as 

R 

Lqeom = ~Tp, j (27) 
167T 

where R is the scalar curvature, but we realize that a different Lagrangian 
might eventually be shown to be more appropriate in a correct theory of 
quantum gravity. 



6 Perfect fluid models 

For a perfect fluid, the energy momentum tensor is 

r = (/3+p)«V + pf, (28) 

where p is the pressure, p is the density, and u is the 4-velocity. For solutions 
to Einstein's field equations for a perfect fluid, (E7J) becomes 



Lgeom = ~P - ~P, (29) 



and we can take the Lagrangian for the matter as [[Schutz and Sorkin 197711 



Lmatter — P (30) 

Thus, the classical action for perfect fluids is 

hi = 1/ (-S (4) ) 1/2 (P - P)d A x ur-V Y U 'i^f<l : '-r. (31) 
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We can take 

V = (7 - 1)P (32) 
for the equation of state, where 1 < 7 < 2. Then (RTT) becomes 



U = \J (-S (4) ) 1/2 (2 - 7W 4 * - ^ / (g (3) ) 1/2 9 m ^d s x. (33) 



Equation (|33"D diverges for a spatially open universe. The significance of 
that might be that only spatially closed universes make sense. On the other 
hand, it might be that the calculation of the action for the correct theory of 
quantum gravity will give a finite value for the action, even for a spatially 
open universe, but here, we shall restrict our calculation to the case of a 
spatially closed universe. 



7 Spatially homogeneous spacetimes 

The integration in (^8j) is an integration over functions g% and 4>i defined on 
Si. In that sense, it is similar to a path integral. For example, there are 
six independent functions that define g%. As in the integration for a path 
integral, there are approximations that can be made to reduce the number 
of integrations that must be performed. 

Here, we want to consider matter distributions similar to that observed, 
at least for the large scale in our universe. Thus, we want to restrict the 
integration in fll8|) to classical spatially homogeneous 4-geometries that have 
a homogeneous matter distribution in calculating the classical action in the 
exponential. The integration in flTjD would then be over the 3-geometries 
that form the boundary of those 4-geometries on Si. 

As an example, we shall use Einstein's General Relativity for the classical 
4-geometries, but the same calculations could be done for other classical 
gravitational theories, in case it turns out that General Relativity is not 
the correct theory of gravity. Thus, we want to consider the integration in 
C P~8[ ) in which the classical 4-geometries used to calculate the action in the 
exponential are restricted to Bianchi cosmologies. 

The appropriate calculation would be to consider the most general Bianchi 
model, with all of the parameters that describe that model, and carry out 
the integration over all of those parameters. We notice that the Bianchi 
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parameters (which are time independent) define the initial three geometry, 
and therefore are valid integration variables in (0). On the other hand, 
if it is suspected that the saddlepoint for the integration will correspond 
to the Friedmann-Robertson- Walker (FRW) model, then one can restrict 
consideration to only those Bianchi models that include the FRW model as a 
special case, and consider integration in (|i~8|) for only one Bianchi parameter 
at a time, holding the others fixed at the FRW value. Here, we do that for 
only one of the Bianchi models for illustration. 

In choosing which Bianchi model to use, we would like one that has 
a parameter that can be varied continuously to give the FRW model. In 
addition, we would like to choose a parameter that represents rotation of 
inertial frames relative to the matter distribution. In that way, we could 
directly test the ability of quantum selection to implement Mach's ideas 
about inertia. 

So far, I have not been able to find a completely satisfactory example. 
Although the Bianchi IX cosmology is often used to represent anisotropy, it 
seems inappropriate for the present case because it is only a superposition of 
gravitational waves on a Friedman-Robertson- Walker background. 

The Bianchi Vlh model seems to be a better homogeneous model that 
has a parameter that represents an angular velocity of inertial frames relative 
to matter, and setting that parameter to zero seems to give the FRW metric. 
However, there seem to be some difficulties with the Bianchi Vlh model being 
able to change continuously into the FRW model, and also a possible problem 
with the topology. Until I find a better example, however, I shall use this 
one. 

We use the solution for the Bianchi Vlh model from [[Ellis and MacCallum 1969 



with h = —1/9. This cosmological model is relevant here because it has a 
relative rotation of inertial frames with respect to the matter. Specifically, 

n(t) = - b - - (34) 

V ; Y 2 {t)Z{t) V ; 

is the angular velocity in the rest frame of an observer moving with the fluid, 
of a set of Fermi-propagated axes with respect to a particular inertial triad. 
The parameter h is an arbitrary constant of the model, and is zero if and 
only if there is no rotation of inertial frames relative to matter. Thus, we are 
interested to know the dependence of the classical action on b. 
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After some algebra, we have 

3vr 2 r*Y(t)Z(t) 



4a J to X(t) 



dt, (35) 



where the spatial part of the 4-volume integration has already been carried 
out, a is a parameter of the model, t corresponds to the surface Si in (|i~6|) 
and (|l8|) and is enough larger than the Planck time T* that the semiclassical 
approximation is valid, the upper limit in ( |35"D corresponds to the surface S2 in 
(|16D and ([Ts|) , and X(t), Y(t), and Z(t) are functions of the model that must 
be determined by differential equations given by [[Ellis and MacCallum 1969 
As expected, the surface term in (|33|) has canceled. 
If we define 



-3k 



3/2 



r"(t) = X(t)Y(t)Z(t) -^— 2 (36) 



and 



1 + "W= U x(jt) w (37) 



then the classical action in ([35|) becomes 

„2\ 1/2 



3tt 2 Z + q l/alY" r (1 + a)r 



-3k 



rii±H£*, (38) 

Jr r 



where f = dr/dt, tq is enough larger than the Planck length L* that quantum 
effects can be neglected, and k = +1 for a closed universe. 
For the h = —1/9 case, we have 

g = -3a (39) 



if and only if b 7^ 0. Substituting ( p9|) into (|38|) gives 

3tt 2 /-1\ 1/2 r (l + a)r 



id 



dr, (40) 



The form of the equation of state in (|32|) allows one of the differential 
equations for the model to be integrated in closed form to give 

8vrp = 3r^- 2 r" 37 (41) 
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where r m is a constant of integration that depends on the amount of matter 
in the universe and the speed of expansion relative to the gravitational at- 
traction. Equation ( [Hp shows that r m is a measure of the amount of matter 
in the universe for a given value of r. Therefore, we might expect Machian 
effects (inertial induction) to increase for larger values of r m . 
Using PH), we have 

. 2 /r m \ 37 ~ 2 , k b 2 r 2 ( a 



f 2 =(^L) -k-ka- '" , —— + — — • (42) 
\rj (2a ) 6 3(l + a) 2 r 4 12U + a/ V ; 

For the isotropic case, only the first two terms on the right hand side of (ff^) 
are nonzero. r m is the value of r where those two terms are equal. For a 
closed universe for the isotropic case, r m is the maximum value of r. 
Equation fl42|) can be written 



r 



2-3 7 k b 2 V 2 



\(r~) - fc(1 + a) - (2a o)6 3(l + a )V + ^ (43) 

and the remaining differential equations to solve are 

k 2b 2 

V ' = 3 * (1 + Q ') r -(2^(TT^ (44) 

and 

—a = ^ (45) 
where ((45|) is the definition of V(t), and r(t) represents the size of the uni- 
verse. The angular rotation ( |34| ) of an inertial frame relative to local matter 
is given by 

/ -k\ 3/2 b 

n W = (^) [TT^Fw (46) 

It is not possible to solve the differential equations exactly in closed form, 
but we can find approximate solutions. In the limit as b approaches zero, 
it is valid to neglect all but the first term under the radical in (f43|). The 
appendix then gives 

1 + a « (1 + a )exp { - 2r ° - - " (rl^ 3 - r| 7 " 3 )r^ §7 + 

I 37 - 6 1 + a 
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(fo - 6)r 3 ^ 2 - (67 - A)rf +1 rh~3 + (3 7 + 2 ) r 3 l 7 ~ 2 
(3 7 + 2)(37-2)(37-6)r^- 2 

8 k b 2 r 3 ^ 6 - 2rf 7 ~ 3 r^~ 3 + rg 7 ~ 6 
~T(2a ) 6 ( 7 -2)(3 7 -6)r^- 2 



12Jfc 



(l + a„){l+ *° (r h-3 _ r b-^b 

{ 37 - 6 1 + a 

(3 7 - 6)r 3 ^ 2 - (67 - 4)rf 7+1 rf 7 " 3 + (3 7 + 2)r 3 1 7 ~ 2 



(37 + 2)(37-2)(3 7 -6)r^- 2 

6 o^i 7_3 ^#7-3 , ^3 7 -6 ' 



8k b 2 r 3 ^ 6 - 2r| 7 ri 7 " 3 + r 3 
~T(2a ) 6 ( 7 -2)(37-6)r^- 2 

(l + a ){l + 



12k 



r 3 7 -2 



(3 7 + 2)(3 7 -2)r^ 2 

8fc b 2 rg 7 ~ 6 I 

3 (2a ) 6 (7-2)(3 7 -6)r 37 - 2 J ' 



(47) 



This solution is valid for r smaller than r m if 6 is small enough. 

We assume that 7 changes from 71 to 72 at r — r%. To satisfy continuity 
of p at r = 7*1, we must have r m change from r m i to r m 2 at r = ri, where 



( 37! -372 
n \ 3 71 -2 

Equation (|4"7| ) applies for r <T\ with 7 = 71 and r m = r ml . 
For r > ri, (flTf) gives 



(48) 



l + a « (1 + ai)ea:p { ^ Ql (r^ 2 3 - rf 2 3 )r^ 2 272 + 

[ 372 - 6 1 + Qfi 

12fe (372 - 6)r 3 ^' 2 - (672 - A)4 12+ \h^ + (3 72 + 2)r^- 2 
(372 + 2)(372 - 2)(37 2 - 6)r^r 2 
8A; 6 2 r 372 " 6 - 2rf 72 " 3 r§ 72 ~ 3 + r? 72 " 6 
~T(2^ (72 - 2)(3 7 2 - 6)rir 2 
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(1+ai) 1 + 



2r 3 



372 - 6 1 + ai 



- n 



!72-3 v l-| T 2 
)'m2 



+ 



12k 



(37 2 - 6)r 3 ^- 2 - (672 - 4)r 1 572+1 ri 72 ~ 3 + (372 + 2)r{ 



372-2 



(372 + 2)(37 2 - 2)(372 - 6)r 



372-2 

Tfl2 



8 A; 6 2 r 372 " 6 - 2r x 272 ?J r b 2 ~ 3 + rf 
_ T(2^ (72 - 2)(3 72 - 6)r^- 2 
« (1 + ai), 

where continuity of a and d at r = ri requires that 

d 



,372-6 



(49) 



1 + «i ~ (1 + a )exp 



2r 3 



771-3 



| 7l -3, l-| 7 i 



371 - 6 1 + a 

,371-2 



ml 



+ 



(3 7 i - 6)rr i_z - (671 - 4)rl 7l+ V 1 f71 - 3 + (3 7 i + 2)r a 



12k 



,371-2 



(3 7l + 2)(3 7 i-2)(3 7 i-6)C 1 1 



371-2 



8/c b 2 r 



371-6 
1 



§71-3 §71-3 

2r rf 



■ 371—6 

+ >0 



3 (2a )6 ( 7l _ 2 )(3 7 i-6)r 3 :r 2 



(l + a )U + 



2r o "0 / §71-3 § 7 i-3 N 1-I71 



12fc 



371 - 6 1 + a 
(37i-6)r 37l - 2 -(67i-4)r 2 '^V 2 



3 7i+l 471-3 



) r ml 

+ (3 7 i + 2)r, 



371-2 




(3 7 i + 2)(37i-2)(37i-6)r^- 2 



8fc b 2 r 



371-6 
1 



§7i-3 §7i-3 

2r rf 



■ 371—6 

+ ?0 



3 (2a ) 6 ( 7l - 2 )(3 7 i-6)r 3 :r 2 
(l + a ){l+ 



12fc- 



r 



371-2 
1 



(3 7 i + 2) (371-2)^ 
8A; 6 2 



371-2 



371—6 



3 (2ao) 6 (7i-2)(3 7 i-6)r^r 2 J ' 
We then perform the integration in (f40"l). The result is 



(50) 



h 



3tt 5 



{-k)V 2 (^) 2 (l + «o) 



37i + 2 



Vr m i/ 



571+1 



(-) 



I71+1 
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24k 



\r m l 



71-1 



/r m iV 

(3 7l + 2)(3 7l -2)(9 7l -2) vr m J 
8A; 6 2 r"t 



§71+1 



r m i 



§71+1 



3 (2a ) 6 (71 - 2) 2 3(3 7 i + 2) 



_ro_ 

'ml 



+ 



1 + 12fc 



r 



371-2^2-371 



ml 



«{. /,2 r -4 r J7i-o r o- 
u ' ml ' 'ml 



371 ~ 6 /rmlV 

Vr m2 / 

371— 6 6— 371 



(3 7a + 2) (3 7l - 2) 3 (2a ) 6 (71 - 2) (3 7l - 6) 



3 72 + 2 



/ r \ §72+1 / ri \ §72+1 

\r m2 J Vr m2 / 



Equation (|5l|) neglects all but the first term under the radical in ([43]). In 
making the calculation, I actually included the other terms under the radical 
to first order, but then determined after the integration that they could be 
neglected. 

Neglecting some small terms gives 



k 

h 



3tt 2 



2 , 2 

(1 + OLQ 



f— y 

3 72 + 2 \r m2 J 



J72+1 



oh, ^2„-4 371-6 6-37J 
_ Ofc 'ml '0 'ml 

~' 3 (2a ) 6 (7i-2)(3 7l -6) 
Substituting (ft§D into (|52]) gives 
3tt 2 



(52) 



id 



2 2 
(1 + a N 



2 : <-*) 1/2 {w) 



(— Y 

3 72 + 2 Vr m2 / 

„37i-6 6-371 
'0 r m2 



772 + 1 



3 Vr m2 7 (2a ) 6 (7i-2)(3 7l -6) 
To get a rough estimate, we take 

7i = 4/3 

to represent a relativistic early universe and 

72 = 1 



(53) 



(54) 



(55) 
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to represent a matter-dominated late universe. Substituting ( |54| ) and (|55|) 
into ([5l|) gives 

37T 2 



Id 

h 



2 

+ 5 



Vr m2 / 

n 

r 

2kb 2 



[1 + Oo) 



/ r \5 /ri \ 5 / 2 / r \ 



r 

'"ml 



Vr m2 / 



r 

2 



I (-) / 

r 



r 

''ml 



^0 

''ml 
2 



/r m iV 



/r m iV 



{2a Q fr 2 ml rl 
3 \r m2 ) 



2 / r \ 5 / 2 / 
5 \r m2 J 1 



r 

r m i. 



_Tp_ 

''ml 



7 



(56) 



Neglecting some small terms, letting k = +1 for a closed universe, and using 
gives 



h 



3i7T 2 /r m2 \ 2 
2 V L* / 



1 + «o) 



1 + 



n 

r 
2b 2 



r 
r 



3 

n 



(-) 



5/2 



J" m2 



+ 



1/2 



f \ 7 / 2 
?"m2 



r \ 



.3/2 



5 (2a ) 6 r m2 rg VrV \r m2 / 

1 + a ) 

2 26 2 



3/2' 



3m 2 /r m2 \ 2 



2 \ L* J 
'2 f r 

5 V?"m2 



5 (2a ) 6 r 2 a2 rg \n/ Vr m2/ / 



3/2' 



(57) 



Because the parameter b is an initial value for the cosmology, it is one of 
the variables of integration in (|H|). In making the saddlepoint approximation 
for that integration, we need to locate the saddlepoint (that is, the value of 
b that makes the action in (|52"1), (|53"D , or (|57|) stationary. We see that the 
action is stationary with respect to variation of b at the isotropic case of 
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b = 0, as expected. The range of values of b that contribute significantly to 
the integral in fllTf ) is given by (0). That is 



Ici{b) I d (b = 0) 



h h 

Thus, substituting into (|5SD gives 



< 1. 



(58) 



3tt 2 /^n2\ 2 (l + a )6 2 / r 
2 U*J (2a ) 6 r^ 2 lr m2 



§72 - 371+1 



16 



( 7 i- 2) 2 9(3 72 + 2) \rj 



r\ 3(71-72) / rm2 \ 



6-371 



The approximations made so far are valid whenever 
Substituting ( 



and (p5) into (|5"9|) gives 



6vr 2 



(1 + ao) 



r \ / r \ f /L* 



rj \r m2 J \r Q J (2a ) 6 L* 4 



< 1. 
holds. 

< 1. 



(59) 



(60) 



This gives 



b< 



V5(2a ) 3 fr^ 1 / 2 



L* 2 /TrraX 

(1 + ao) 1 / 2 V r J 



7T\/6 

Thus, from (46|), the rotation rate of inertial frames is 



L* 



(61) 



\m\ < 



i 



VE f L* \ 2 

^/Q\r~J (l + a ) 1/2 [l + «(t)]r m2 \r{t) _ 



r m 2 



IS 

4 



r J 



1/2 



L* 



(62) 

If we now take the Planck length L* to be 1.6 x 10 33 cm, use the Hubble 
distance of 1.7 x 10 28 cm for r m2 , a tenth of that for r, neglect a and «o 
compared to 1, and take 

r ! = TKK ( 63 ) 



as an estimate that the universe changed from radiation-dominated to matter- 
dominated when the universe was about one-hundredth of its present size 



Weinberg 1972| , Section 15.3, p. 481], then we get 

|0[ < 1.6 x 10~ 130 radians per year, 



(64) 
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which is much less than the bound set by experiment of 10~ 14 to 7 x 10~ 17 
radians per year if the universe is spatially closed [ Hawking 196S . 

The rotation rate in ( |5^ ) and QS1]) is so small because the Planck length 
is so much smaller than the Hubble distance. 

That the small value of allowed rotation rate depends mostly on the 
universe being much larger than a Planck length rather than on details of 
the model suggests that the result has some generality. 

We notice also, that the selection criterion in ( pljj) is so sharp that the 
initial wave function in the integration in fll8|) would have to be very sharply 
peaked to overcome it. 



8 Discussion 

We see that considerations of quantum cosmology show how a range of clas- 
sical cosmologies can be selected that contribute significantly to the wave 
function in the final state. The effect enters through the action. Using 
semiclassical calculations gives results that should not depend on particular 
details of the theory of quantum gravity. 

For our universe (which is much larger than the Planck length) the selec- 
tion is very sharp. The initial wave function over 3-geometries would have to 
be extremely sharp (not a probable occurrence) to dominate over the effect 
of the action. 

The selection process seems to occur very soon in the development of 
a cosmology. That is, for a broad wave function over 3-geometries in the 
initial state, the wave function becomes sharply peaked after the universe 
has become a few orders of magnitude larger than the Planck length. 

A different choice than (|3^) [[5 chut z 1976| is 



■'matter 



p. (65) 



The choice in (|65| ) gives a third of ([H]) for the total action. A correct theory 
of quantum gravity will determine which (if either) of these two choices is 
correct, but for this illustration, a factor of three in the action makes little 
difference. 

It appears likely now, however, that there is not enough matter to keep 
our universe from expanding forever (e.g. [|Colcs and Ellis 1994|| ). To accom- 



modate that with a spatially closed universe within General Relativity would 
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require a positive cosmological constant (e.g. [[Hawking 1998|| ). I shall try to 
include the cosmological constant in a future calculation. 
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A Bianchi VI h Models 

We start with Equations (6.7) of Ellis and MacCallum (1969). For the case 
of zero cosmological constant, A, these can be written as 



P 



P + P 



3R 



4vr(p -p) 
Att(p-p) = 



R 



R- 



-3 



XX 2{a 2 + q 2 ) 



+ 



2b 2 



4 



X 2 Y±Z 2 
2{a 2 + a q ) 2b 2 



X 2 



Y A Z 2 



8vrp + 
where 



Att{ P -p) = R- 3 [R 3 

3ag + ql = 1 
X 2 2 





2b 2 



Y 4 Z 2 ' 



R(t) 3 = X(t)Y(t)Z(t), 



(66) 

(67) 
(68) 
(69) 
(70) 

(71) 



b, ao, and go are constants that are parameters of the model, and the 
variation of p, X, Y, and Z with time is determined by fl66"D through ([701). 

According to Ellis and MacCallum (1969), (|70|) is a first integral of the 
other equations. If I understand that correctly, then taking the derivative of 
(|70|) should be a combination of the other equations. When I try that, the 
result is close, but not quite correct. I have not been able to find a similar 
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equation that does work, although I have found one that works for the special 
case (go + 3ao = and XY = Z 2 ) that I use later. This is 



8vrp + 



3ag + ql = 1 
X 2 2 



± 



Y_ 



(72) 



We can add 
(A+B+C) 
where 



8|), and ( |69|) with coefficients A, £>, and C to give 



4^r(p-p) + ^f 



(73) 



We can take special cases of A, B, and C. Taking A = B = C = l/3in 
gives 

A,{p- P ) = R-\R 2 ky-f^±^, 

where I have used ([71]). For another special case, we take A = B = 1 and 
C = -2 in (fF3) to give 



(75) 



3a + q 
- 2 go x2 ■ 



where 



U 



XY 

'z 2 ' 



(76) 



(77) 



For a third special case, we take A = —4/3 and £> = C = 2/3in (|73|) to give 



iT 3 R 3 



a 



Ab 2 8q 2 



where 



1 + aJ Y 4 Z 2 3X 2 ' 

y2/3 Z 2/3 



1 + a 



(78) 



(79) 



X 4 /3 • 

We can use (|TTD, ([77D, and QT9"D to determine X, Y, and Z in terms of R, U, 
and a. This gives 



R 

(1 + a) 1 / 2 ' 



(80) 
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and 



Y = U 1 ' 3 (l+a) 1 ' 2 R, 

z- R 



Using ©,(13), and © in and © gives 



IT- 



-2?, 



iT 3 i? 3 



a 



3R 2 

(3a + g )(l + a) 
46 2 8g 2 (l + a) 



l + a7 (1 + a) 2 U 2 / 3 R e 



and 



i2. 



1 a 



21 + a 



+ -x 



1 a 



3R 2 



np + (3a 2 + q 2 )(l + a) 



21 + a/ \3UJ \3UJ 
b 2 



3 ' 3R? 
If we change variables from R to r, where 



3{1 + a) 2 U 2 / 3 R e ' 



-3k 



2 R, 



V 3«o + % 

(and where k is +1 for a closed universe and -1 for an open universe) 
(|83"D through (j86l) become 



r (r r)' = 47r(p — p) 



2fc(l + a) 



-3 / ^3 



f/Y Qkq (3a + q )(l + a) 



U 



(3a 2 + q 2 )r 2 
Ab 2 



"3 „3 



a 



1 + a 



-3k 



+ 



8kq 2 (l + a) 



3al + ql) (1 + a) 2 U 2 / 3 r 6 (3a 2 , + %V 
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and 



1 d 



21 + a. 
np k(l + a) 



1 a \ aca fur 

21 + a) \3U) + \3U 



+ 



-3k 



3a 2 + q 2 ) 3(l + a) 2 f/ 2 / 3 r 6 ' 



If we define W, V, and K by 



and 



r 2 r = W, 
a 2V 



1 + a r- 
U _3K 



then (^8|) through ([H]) become 

W = 4vr(p - p)r 3 - 2fc(l + a)r, 
2kq (3a + q )(l + a)r 



V 



K 



-3k 



3a 2 , + 
2b 2 



+ 



4kq 2 (l + a)r 



3al + ql) (l + a) 2 f/ 2 / 3 r 3 3a§ + g 2 ' 



V 2 + VK + K 2 87ipr 6 



W 

3 

-k(l + a)r 4 + 

If we take 

for the equation of state, where 



-3k 



, 3 b 2 



3al + ql) 3(l + a) 2 U 2 / 3 ' 
P = (7 - 1)P 



1 < 7 < 2 
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is a constant, then can be integrated to give 

8np = 3r^- 2 r" 37 



(101) 



where r m is a constant of integration. Substituting (p9|) and (|101|) into (|95]) 
and gives 



= 3(1 - 7/2)r^ 7 - 2 r 3 - 37 - 2k(l + a)r, 



(102) 



V 2 + VK + if 2 



-&(l + a)r 4 + 



_i_ „37-2 6-37 

' 771 



-3Jfe 



Zal + ql) 3(1 + a) 2 f/ 2 / 3 ' 
Let us now convert to a set of dimensionless variables defined by 

r = t/r m 
x = r/r m 

y = v/rl 

z = W/rl 
C = K/r 2 m 
D = b/r 2 m 

Let us also use ' for d/dr. Then © through (0), ([Top , and ( p3|) 



x 2 x' 



a 



1 + a 



2y 
,r 3 



3C 



C" 
-3A; 



2A;g (3ao + <?o)(l + a)^ 



3a 2 , + <ll 
2D 2 



+ 



4/cg 2 (l + a)x 



3a 2 + qiJ (l + a) 2 U 2 / 3 x 3 3a 2 + ql 



(103) 

(104) 

(105) 
(106) 
(107) 
(108) 
(109) 
become 

(110) 
(111) 
(112) 
(113) 

(114) 
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z ' = 3(1 - 7 /2)x 3 ~ 37 - 2fc(l + a)x, (115) 

For the h = —1/9 case, we have 

go = -3a . (117) 

if and only if b ^ 0. However, the 6 = case is only a single point. When 
integrating over b, the behavior for small b is more important than exactly 
at b = 0. Therefore, we shall use ( 117 ) in any case. Therefore, from ( |113| ), 



we have that C is a constant. 

So both b and C are constants that are determined by initial conditions. 
There are two possibilities. Either they are connected by a relation, or they 
are independent. There seems to be no obvious connection between them 
from the equations, so we shall assume they are independent unless we shall 
find it necessary to make a connection to get a consistent solution to the 
equations. Therefore, we can assume that both b and C must be integrated 
over on the initial hypersurface. For this demonstration, however, it is suffi- 
cient to integrate over only one initial constant, b. Therefore, we shall fix C 
at the value we would guess for the FRW case, which is zero. Therefore, we 
take C to be zero. 

That means from (|112[) that U is constant. Again, we guess that U is 
independent from b, and choose the isotropic value. In the isotropic case, 
we would have X = Y = Z, and therefore, from ( |77| ) we have U = 1 as the 
isotropic value. 



Therefore, from (114) we have 



k 2D 2 

y' = m + <x)x- {2aoril + a)2x3 , (us) 



and from ( 116 ) we have 



;W--- M l + ^-(^3U^4. (U9) 

We notice at this point that if we take the derivative of ( |119| ) and substitute 
from ( |110|) , (|1 1 1|) , (p.l8|) , and (|115 ), that we get an identity, confirming the 



consistency of the equations at this point. 
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Combining ( |110| ) with ( |119| ) gives 



x 



\ 



x 



2-37 



k(l + a) 



(2a ) 6 3(l + a) 2 x 4 3x 4 



+ (120) 



Combining ( 118 ) with ( P-20| ) gives 

3k(l + a)x- k 



dy_ 

dx 



2D 2 



(2a ) 6 (l+a) 2 x 3 



X 



2-37 



k(l + a) 



k D 2 , y 2 

(2a ) 6 3(l+«) 2 a; 4 ~ t ~ 3x 4 



(121) 



Both b and do are constants that are determined by initial conditions on 
the initial hypersurface. They are either related or independent. We assume 
first that they are independent. In that case, the third term under the radical 
in ( p. 20] ) and ( |121|) will get smaller as D gets smaller. We shall assume that 
we can neglect that term relative to the first term under the radical for all 
values of x. We can check that assumption later. There is also the possibility 
of iterating later by assuming that this term is small instead of zero. 

We shall also neglect the last term under the radical in ( |120| ) and ( |121| ). 
We can check that approximation later. There is also the possibility of iter- 
ating by substituting an approximate solution for y into (|121|). We shall see 



that it it permissible to neglect the fourth term for all values of x within the 
integration range if b is small enough. 

Neglecting the third and fourth terms under the radical in ( 120|) and ( |121| ) 



gives 



and 



x 



x 



2-37 



fc(l + a) 



2D A 



dy_ _ 3/c(l + a)x - (2qq)6 (1+a)2x 3 
dx 



-k(l + a) 



We shall assume that 



a < 1. 



(122) 



(123) 



(124) 



We can test that approximation from the solution later and iterate if neces- 
sary. In that case, ( |120| ) and (|121|) become 



X 1 = yV-37 _ 



(125) 
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and 



dy 
dx 



3 kcc 



2D 2 



(2a ) 6 



(126) 



^2-37 _ k 

We notice that ( |1 26| ) could be integrated numerically to obtain an ap- 
proximate solution for y(x). For a closed universe (which is the case we 
are considering), we can get an estimate of that solution, at least for small 
x. The second term under the radical in (|125|) and ( |126| ) is smaller than the 
first term under the radical except at the point of maximum expansion, when 
they are equal (when x is one). Thus, before the universe gets too close to 
the point of maximum expansion, we can neglect the second term under the 



radical in ( |125| ) and ( |126|) to give 



x 



x 



l-§7 



and 



where 



^ = Zkxh - 2kE 2 xh~\ 
dx 



E = D/(2a ) 



(2a ) 3 



We can integrate ( |128| ) to give 



V = 2/o + 



6k 



3 7 + 2 



3 

(X 2 



7+1 rr2~1 +1 - 
Xn 



4k E 2 

y~2' 



[X 2 



§7-3 



Xr 



7-3 N 



(127) 
(128) 
(129) 

(130) 



where, without loss of generality, the constant of integration has been chosen 
such that y = yo when x = x$. 

We notice that on substituting ( |130|) into ( |120| ) and ( |121| ) the terms it 
gives that do not involve h have a higher power of x than the dominant term. 
Therefore, for small b and small x, it was justified to neglect the fourth 
term under the radical in (|120|) and (|121|) , since the terms that involve b are 
proportional to a positive power of b. 

Using (|130|) in ( |1 1 1| ) gives 



a' 



1 + a 



2y x-' 3 + 



12A; 
3 7 + 2' 



[x 2 



-2 

Xn 



§7+1 



X 



8k E 2 

i~2' 



(X 2 



7—6 



— Xr 



7-3. 



X 



(131) 
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Combining ( |131| ) with ( |127| ) gives 

dln(l + a) _ 2 3 7 - 4 , 12fc f 37 _ 3 | 7 +i a 4 8A; E 2 37 _ 7 §7-3 3 7 _ 4 

(132) 

We can integrate ( |132j ) to give 

1 + a 37 — 

12k (3 7 - 6)ar^- 2 - (67 - 4)x| 7+1 xi 7 " 3 + (3 7 + 2)x^ 7 " 2 



3 7 + 2 (3 7 - 2)(3 7 -6) 

8k E 2 x 3 ^ 6 - 2xp~\h-z + x 3 ^- 6 



(133) 



3 7-2 37-6 

where, without loss of generality, the constant of integration has been chosen 
such that a = a when x = x . 

Writing (PPQ and (|I33D in original variables using QTUoD , (|TUB| ), ([Tgg), 
and (|9~3] ) gives 

^ = r 3 a _ r 3 d | 
21 + a" 21 + a 

[bk 3 7 + 2 T(2^F 7^2 Jrm ' (134) 

l + a ~ (1 + ao )exp { - 2r ° 3 - - d ° (r§ 7 ~ 3 - rl 7 ~ 3 )r^" |7 + 



12fc 



37 — 6 1 + a 

(3 7 - 6)r 3 ^ 2 - (67 - A)rp +1 rh- 3 + (3 7 + 2)r 37 ~ 2 
(37 + 2)(3 7 -2)(37-6)r 3 ra 7 " 2 

8 k b 2 r 3 ^" 6 - 2r| 7 " 3 rf^ 3 + r 3 , 7 ' 6 
"3"(2a ) 6 ( 7 -2)(37-6)r 37 " 2 



(135) 
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